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We calculate the critical value of the hopping parameter, K^, in Lattice QCD, up to two loops in 

perturbation theory. We employ the Sheikholeslami-Wohlert (clover) improved action for Wilson 

fermions and the Symanzik improved gluon action for 4- and 6-link loops. 

The quantity which we study is a typical case of a vacuum expectation value resulting in an 

additive renormalization; as such, it is characterized by a power (Unear) divergence in the lattice 

spacing, and its calculation lies at the limits of applicability of perturbation theory. 

Our results are polynomial in csw (clover parameter) and cover a wide range of values for the 

Symanzik coefficients c, . Furthermore, the dependence on the number of colors and the number 

of fermionic flavors Nf is shown explicitly. In order to compare our results to non perturbative 

evaluations of fQ coming from Monte Carlo simulations, we employ an improved perturbation 

theory method applied to improved actions. 
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1. Introduction 

In the present work, we study clover fermions [[T]] and Symanzik improved gluons [Q], for the 
calculation of the critical value of the hopping parameter, jq. in perturbative Lattice QCD, to two 
loops order in a perturbative expansion. 

Concerning gauge fields, we employ the Symanzik improved action The purpose of in- 
troducing such an improvement of the lattice action, is to minimize the consequences of a non-zero 
lattice spacing, order by order in perturbation theory. Several choices are made for the coefficients 
of the action corresponding to values which are most often used in the literature. 

The discretization of the theory via lattice regularization introduces some difficulties that do 
not exist in the continuum theory. To be able to recover the continuum hmit, we must demand strict 
locality and absence of doublers. It is known that these requirements lead to breaking of chiral 
symmetry. Of course, while approaching the continuum limit we expect to recover chirality. This 
is the point where one must introduce the hopping parameter. The main idea is to ensure chiral 
symmetry, by setting the renormalized fermionic mass (ms) equal to zero. Because of additive 
renormalization, setting the bare fermionic mass equal to zero is not enough. Hence, there is a 
critical mass, the role of which is to guarantee that tur vanishes. This quantity is directly related to 
the hopping parameter K. Its critical value {Kc) is responsible for restoring chiral symmetry. 

Previous works on the hopping parameter and its critical value appear in the literature for 
Wilson fermions-plaquette action gluons [^] and for clover fermions-plaquette action gluons [Q]. 
The procedure and notation in our work is the same as in the above references. 

Our results for Kc (and thus for the critical mass) depend on the number of colors {N) and 
fermionic flavors {Nf). Besides that, there is an explicit dependence on the clover parameter csw 
which is kept as a free parameter. The dependence of the results on the choice of Symanzik coeffi- 
cients cannot be given in closed form; instead, we present it in a list of Tables and Figures. 

The rest of the paper is organized as follows: In Sec. |2| we define the actions in their discretized 
form, as well as the connection between the hopping parameter and fermionic mass. Furthermore, 
there is a description of our calculations along with the necessary Feynman diagrams. In Sec. |3|our 
results are presented and compared with previous Monte Carlo simulations. Finally, in Sec. ^ we 
apply to our 1- and 2-loop results an improvement method, proposed by us [^, ^, ^J. This method 
resums a certain infinite class of subdiagrams, to all orders in perturbation theory, leading to im- 
proved results. A full write-up of this work, including detailed tables of results is forthcoming [^. 



2. Formulation of the Problem 



We begin with the Wilson formulation of the QCD action on the lattice, with the addition of 
the clover (SW) [|l|] term for fermions. In standard notation, it reads: 

Sl = 5G + ££(4r + m)i|ir/(x)i/A/(x) 

/ ^ 



EI 



+ T^swE E Vfix)(y^vF^vix)Yfix) 

f x,fi.v 



(2.1) 
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1 ^ 

where: F^v = {Q^v - Qvp) (2.2) 

and: Qjxv = Ux,x+jiUx+^.x+^i+vUx+^+v,x+vUx+v,x~\~Ux,x+vUx+v,x+v-^Ux+v-^.x-^i^x-ji,x 

+ Ux^x-ixUx-ji,x-ii-vUx-ji-v,x-vUx-v,x + Ux^x-vUx-v,x-v+iiUx-v+ii,x+iiUx+jx,x (2.3) 

The clover coefficient csw is treated here as a free parameter; r is the Wilson parameter; / is a 
flavor index; a^y = {i/2)[y^, jy]. Powers of the lattice spacing a have been omitted. 

Regarding gluons, we use the Symanzik action, involving Wilson loops with 4 and 6 links: 



2 



CO £ Re Tr ( 1 - C/piaq) + ci ^ Re Tr ( 1 - f/,ect) 

plaq rect 
-C2 £ ReTr (1 - U^^^) + C3 ^ ReTr (1 - U^^^i) 

chair paral 



(2.4) 



The correct classical continuum limit requires: cq + 8ci + 16c2 + 8C3 = 1. The full action is: 
S = Sl + Sgf + 5^;, + Sm; where 5^/, Sgh, are standard gauge fixing, ghost and measure terms. 

The bare fermionic mass mg must be set to zero for chiral invariance in the classical continuum 
Umit. The value of the parameter csw and of the Symanzik coefficients c; can be chosen arbitrarily; 
they are normally tuned in a way as to minimize t&{a) effects. Terms proportional to r in the action, 
as well as the clover terms, break chiral invariance. They vanish in the classical continuum limit; at 
the quantum level, they induce nonvanishing, flavor-independent corrections to the fermion masses. 
Numerical simulation algorithms usually employ the hopping parameter, 

^ = 7^ (2.5) 

Imsa + or 

as an adjustable quantity. Its critical value, at which chiral symmetry is restored, is thus l/8r 
classically, but gets shifted by quantum effects. 

We denote by dm the perturbative contribution that must be added to the bare mass, in order 
to lead to zero renormaUzed mass. At tree level, mg = 0. 



dm = i/m(i_ioop) + dm(^2-\oop) (2.6) 

Two diagrams contribute to Jm(i_ioop), shown in Figure 1. The quantity <i'«(2-ioop) receives 
contributions from a total of 26 diagrams, shown in Figure 2. Genuine 2-loop diagrams must be 
evaluated at mg 0; in addition, one must include to this order the 1-loop diagram containing an 
ff{g^ ) mass counterterm (diagram 23). Certain sets of diagrams, corresponding to renormalization 
of 1-loop propagators, must be evaluated together in order to obtain an infrared-convergent result: 
These are diagrams 7+8+9+10+11, 12+13, 14+15+16+17+18, 19+20, 21+22+23. 



3. Numerical Results 



We have selected a set of most widely used values for the Symanzik coefficients, shown in 
Table 1. (In all these cases, C2 = 0. In general, for given values of Ci = C2 + cj, C2 = ci — C2 — C3 
the dependence on C2 is polynomial and thus we need not choose a numerical value for it.) 
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Action 




CO 


ci 


C3 


Plaquette 




1.0 


0.0 


0.0 


Symanzik 




1.6666667 


-0.083333 


0.0 


TILW, pco = ! 


8.60 


2.3168064 


-0.151791 


-0.0128098 


TILW, Pco = : 


8.45 


2.3460240 


-0.154846 


-0.0134070 


TILW, jSco = : 


8.30 


2.3869776 


-0.159128 


-0.0142442 


TILW, jSco = : 


8.20 


2.4127840 


-0.161827 


-0.0147710 


TILW, jSco = : 


8.10 


2.4465400 


-0.165353 


-0.0154645 


TILW, /3co = 1 


8.00 


2.4891712 


-0.169805 


-0.0163414 


Iwasaki 




3.648 


-0.331 


0.0 


DBW2 




12.2688 


-1.4086 


0.0 



Table 1: Input parameters cq, ci, C3. 




3 4 5 6 7 8 




9 10 11 12 13 14 




15 16 17 18 19 20 

Figure 1: One-loop diagrams 



contributing to <i»i(i-ioop)- 
Wavy (solid) lines represent 
gluons (fermions). 




Figure 2: Two-loop diagrams contributing to fiwp-ioop)- Wavy 
(solid, dotted) lines represent gluons (fermions, ghosts). Crosses de- 
note vertices stenmiing from the measure part of the action; a soUd 
circle is a fermion mass counterterm. 
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The contiibution of the I 1-loop diagram to dm, can be expressed by: 

di=^^^s'-h's^4^ (3.1) 

where ej'^ are numerical constants and their values depend on the Symanzik coefficients c,-. The 
dependence on csw is seen to be polynomial of degree 2 (/ = 0, 1, 2). 

The contribution of 2-loop diagrams without closed fermion loops takes the form 

^/ = ^^^/-I4wA^^4^P'^ (3.2) 

where the index / runs over all contributing diagrams, j = 0,2 and ^ = 0, 1, 2. The coefficients 
ef'^''^'^ exhibit a further dependence on Symanzik coefficients (only through the combinations Ci , C2), 
which cannot be expressed in closed form and is presented numerically in what follows. 
The contribution of 2-loop diagrams, containing closed fermion loops, has the form 

^/ = ^^^A^//-i4w# (3.3) 

In order to enable cross-checks, numerical per-diagram values of the constants e\'\ and 
ef^ are presented in our forthcoming publication for the Iwasaki action. 

The total contribution of 1-loop diagrams, for N = 3, can be written as a polynomial function 
of the clover parameter csw- For the Wilson and Iwasaki actions we find, respectively: 

^'"(I'-Zp) =S^{- 0.434285489( 1 ) + 0. 1 159547570(3) csw + 0.0482553833 (1 ) c^w) (3.4) 

^"^O-'ioop) (-0.2201449497(1) +0.0761203698(3) csw +0.0262264231(l)ciw) (3.5) 
To illustrate our 2-loop results for some particular choices of the action, we set A'^ = 3, C2 = 
and we use three different values for the flavour number: Nf = 0, 2. Thus, for the Wilson action: 

Nf = 0: t/m(2_ioop) = / ( -0.1255626(2) + 0.0203001 (2) csw +0.00108420(7) c^w 

-0.00116538(2)c|w -0.0000996725(l)4w) (3.6) 

Nf = 2: dm(^2~ioop) 0.1192361(2) +0.0173870(2) csw +0.00836498 (8) c^w 

-0.00485727(3) ciw -0.001 1561947(4) c^w) (3.7) 

and for the Iwasaki action: 

Nf = 0: t/m(2_ioop) = / (-0.0099523(2) - 0.0024304(5) csw - 0.00232855(4) ciw 

-0.00032100(2) ciw - 0.0000419365(1) c|w) (3.8) 
Nf = 2: t/m(2_ioop) = (- 0.0076299(2) - 0.004073 1(5) csw +0.00102758 (6) c|w 

-0.00242924(3) 4^ - 0.000457690(2) c^w) (3.9) 

In Figures 3 and 4 we present the values of dm^2-\oop) for = 0, 2, respectively; the results are 
shown for all our choices of Symanzik actions, as a function of csw (A'^ = 3, C2 = 0). 
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g -0.06 

S -0.08 - 




0.0 0.5 1.0 1.5 2.0 

Figure 3: Total contribution of 2-loop diagrams, 
for N = 3, Nf ^0 and C2 = 0. Actions (top to 
bottom): DBW2, Iwasaki, TILW(j3=8.00, ), 
8.20, 8.30, 8.45, 8.60), Symanzik, Plaquette. 

4. Improved Perturbation Theory 



0.02 
0.00 




Figure 4: Total contribution of 2-loop diagrams, 
for N ^ 3, Nf ~ 2 and C2 = 0. Actions (top to 
bottom): DBW2, Iwasaki, TILW(j3=8.00, 
8.20, 8.30, 8.45, 8.60), Symanzik, Plaquette. 



We now apply our method of improving perturbation theory based on resummation of 
an infinite subset of tadpole diagrams, termed "cactus" diagrams. In Ref . [0] we show how this 
procedure provides a simple way of dressing (to all orders) perturbative results at any given order 
(such as the 1- and 2-loop results of the present calculation). Some alternative ways of improving 
perturbation theory have been proposed in Refs. |^ |l^. In a nutshell, our procedure involves 
replacing the original values of the Symanzik and clover coefficients by improved values, which 
are explicitly computed in |^. Taking also due care to avoid double counting of diagrams, we 
calculate the improved ("dressed") value <im* of the critical mass {N = 3, C2 = 0). 

We choose to study the case of the Wilson action (j3 = 5.29) and the Iwasaki action (j3 = 1.95) 
with N = 3 and Nf = 2. Using these values, the contribution to dm'l^^-ioop] ^ polynomial of csw- 

^'"(2-ioop),wiison = - 0.77398(5) +0.1632999(3)csw + 0.06225(2)ciw 

- 0.00440064(8) c^w- 0.000737797(1) 4w (4-1) 
^'"(2-ioop),iwasaki = - 0.0813302(9) + 0.043030(3) Csw + 0.0308196(2) c^^ 

- 0.00767090(8) c|w- 0.001160923(1) c^w (4-2) 

The comparison between the total dressed contribution <im* = dm'^^ ^^^^^ +'^"^*-ioop) ^^'^ 
unimproved contribution, dm, for the plaquette action (j8 = 5.29, Nf = 2) is exhibited in Figure 5, 
as a function of csw- Similarly, dm for the Iwasaki action (/3 = 1.95, Nf = 2) is shown in Figure 6. 



Action Nf p csw K:i-loop K"2-loop K"i*ioop '^2-loop """^ 

Plaquette 6m 1.479 0.1301 0.1335 0.1362 0.1362 0.1392 

Plaquette 6.00 1.769 0.1275 0.1306 0.1337 0.1332 0.1353 

Plaquette 2 5.29 1.9192 0.1262 0.1307 0.1353 0.1341 0.1373 

Iwasaki 2 1.95 1.53 0.1292 0.1368 0.1388 0.1379 0.1421 

Table 2: 1- and 2-loop results, and non-perturbative estimates for k^- 
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Finally, in Table 2, we present a comparison of dressed and undressed results with non perturbative 



estimates for K^^^ [[111, |12|, [13|]. We observe that improved perturbation theory, applied to 1-loop 
results, already leads to a much better agreement with the non perturbative estimates. 





Figure 5: Improved and unimproved values of 
dm as a function of csw, for the plaquette action. 



Figure 6: Improved and unimproved values of 
dm as a function of csw. for the Iwasaki action. 
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